Percolation in a kinetic opinion exchange model 
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We study the percolation transition of the geometrical clusters in the square lattice LCCC model (a 
kinetic opinion exchange model introduced by Lallouache et al. in Phys. Rev. E 82 056f 12 (2010)) 
with the change in conviction and influencing parameter. The cluster comprises of the adjacent sites 
having an opinion value greater than or equal to a prefixed threshold value of opinion {yi). The 
transition point is different from that obtained for the transition of the order parameter (average 
opinion value) found by Lallouache et al. Although the transition point varies with the change in the 
threshold value of the opinion, the critical exponents for the percolation transition obtained from 
the data collapses of the maximum cluster size, cluster size distribution and Binder cumulant remain 
same. The exponents are also independent of the values of conviction and influencing parameters 
indicating the robustness of this transition. The exponents do not match with that of any other 
known percolation exponents (e.g. the static Ising, dynamic Ising, standard percolation) and thus 
characterizes the LCCC model to belong to a separate universality class. 

PACS numbers: 64.60. ah, 05.50.+q, 05.70.Fh 



I. INTRODUCTION 



Geometrical percolation transition has been a long 
studied subject [HQ- It is characterised by a set of 
universal critical exponents, which describe the fractal 
properties of the percolating medium at large scales and 
sufhciently close to the transition. The exponents only 
depend on the type of percolation model and on the spa- 
tial dimension. The occupancy of the sites or bonds of a 
percolating system is controlled by a parameter and at a 
critical value of that parameter the cluster sizes (defined 
by the number of adjacent sites posessing a pre-defined 
common feature) goes to infinity which we call percola- 
tion transition. This phenomenon has been extensively 
studied for the thermal excitation of the two dimensional 
Ising model and in this case the system undergoes perco- 
lation transition [3|4q| at the same critical temperature as 
the magnetization [3, Q ■ In case of Ising model, the geo- 
metrical cluster is defined by the adjacent sites consisting 
of parallel spins. The transition point differs in case of 
higher dimensions [H, Hj. The percolation exponents of 
the geometrical clusters are identical for the models be- 
longing to the same universality class (Ising and Z(3) 
symmetric models [13, Ell)- Recently the dynamical per- 
colation transition has been studied for 2d Ising model 
by applying pulsed magnetic field [l2|. The critical ex- 
ponents were different from that of the static percolation 
transition associated with the thermal transition of the 
Ising model indicating a different universality class. The 
distinct crossing point of the Binder cumulant of the or- 
der parameter for different system sizes at the transition 
point also characterizes the being of the dynamical peco- 
lation transition in a different universality class. 
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Study of social dynamics has been very popular in re- 
cent times and to capture the basic idea of consensus 
formation concepts of statistical physics has been applied 
largely A large number of models has been stud- 

ied (voter model [l5l.Tl6t. Sznajd model [l^ etc) so far. 
In some models opinions have been considered as a con- 
tinuous variable [18.-^23i]. The spreading of an opinion 
through the society may be compared with the percola- 
tion problem in physics and has been studied for non- 
consensus opinion model earlier j24i] and was found to 
belong to the same universality class as the invasion per- 
colation. In this paper we have studied the percolation 
transition of geometrical clusters in a recently proposed 
opinion model called the LCCC model [13, [l^] in which 
individuals exchange opinions controlled by an influenc- 
ing parameter and a conviction parameter, the values of 
which are equal. The opinion of an individual is taken 
uniformly between —1 and which changes by binary 
interactions, where an individual stays with his own opin- 
ion upto a certain fraction A and takes a random frac- 
tion of a part of another agent's opinion determined by 
the same parameter (detailed discussion has been given 
in the next section). By Monte Carlo simulation it was 
found that below a critical value (Ac « 2/3) of the convic- 
tion parameter the average opinion value remains zero, 
whereas above the critical value the average opinion value 
becomes non-zero. Some critical exponents characteris- 
ing the transition in LCCC model and some variants of 
the LCCC model were studied numerically [IBl- A gener- 
alised version of this model was introduced by Sen [23| in 
which the influencing parameter and the conviction pa- 
rameter were different. A discrete version of the LCCC 
model has also been studied [28| . 

In this paper we have investigated the percolation tran- 
sition of the geometrical clusters of the LCCC model as- 
suming individuals are located on the sites of a square 
lattice. We have defined clusters as a group of adja- 
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cent sites with opinion value equal to or above a preas- 
signed threshold value (ft). The cluster sizes are con- 
trolled by the influencing parameter A and for a fixed Q, 
at a critical value of the influencing parameter A^, the 
percolation transition occurs. We determine the critical 
exponents by finite size scaling analysis of the maximum 
cluster size. The value of the critical point decreases 
with decrease of fl and coincides with that for the tran- 
sition point of Ac = 2/3 (at which the average opinion 
diverges) as f2 — 0.0. But the critical exponents re- 
main unaltered with change in n and also differ with the 
exponents known for the previously known models indi- 
cating the square lattice LCCC to belong to a separate 
university class. We have also investigated this percola- 
tion transition in case of generalised LCCC model (27| . 
where the conviction parameter (A) is different from the 
influencing parameter (/i) and once again found that al- 
though the critical point shifts depending on the values 
of fl, A and /i, the critical exponents remain same. 

The paper has been organised in the following manner 
: In Sec II we give a brief description of both the LCCC 
model and the generalised LCCC model. In Sec III we 
present the description of clusters and measure the crit- 
ical exponents for the square lattice LCCC model. In 
Section IV we measure the exponents for the generalised 
LCCC model and flnally in Sec V we have some discus- 
sions regarding this transition and some conclusions. 



II. A BRIEF DESCRIPTION OF THE LCCC 
AND GENERALISED LCCC MODEL 

The origin of this model is a multi-agent statistical 
model of closed economy [2^ where N agents exchange a 
fixed wealth through pair- wise interaction controlled by a 
"saving" parameter. Lallouache et. al. [H, [25| proposed 
a similar multiagent model to describe the dynamics of 
opinion formation. The basic difference in this model is 
that there is no constraint regarding the conservation of 
opinion. Let there are N agents and each agent i begins 
with an individual opinion Oi G [—1,-1-1]. They exchange 
opinions between each other by binary interactions as 
follows : 
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FIG. 1: (Color online) Maximum cluster size as a function 
of the conviction parameter for four different system sizes 
(L = 60, 80, 100 and 200) for the LCCC model i.e. = A 
and threshold opinion value f2 = 1.0. 

model was as follows [27j : 

o,(t + l) = Ao,(i) +e^Oj(i) 

Oj(t + l) = Aoj(t)+eVKi) (2) 

where A is the conviction parameter and [i, is the influ- 
encing parameter with oi € [—1,-1-1]. The special case of 
A = /J, is the LCCC model. The order parameter is the 
average opinion O = \ X^i^ij/iV. Numerical simulations 
show that the system stabilises into two possible phases 
: for any A < Ac, Oi = Vz, while for A > Ac, O > and 
O — ^ 1 as A -H* 1. In LCCC model Ac « 2/3 is the critical 
point. In case of the generalised model Ac depends on the 
value of and the mean field phase boundary is given by 
A = l — fi/2. If we study these models on a square lattice, 
then also the critical points do not change. Some critical 
exponents characterising the transition in LCCC model 
and some variants of the LCCC model were also studied 
numerically [26] . 

III. PERCOLATION ON SQUARE LATTICE 
LCCC MODEL 



o,{t + l) ^ X[o^{t) + eOj{t)] 

Oj{t + l) = X[oj{t) + e'o,{t)] (1) 

where e, e' are drawn randomly from uniform distribu- 
tions in [0, 1]. In this model the opinions are bounded i.e., 
— 1 < Oi < +1 for all i. Here the parameter A is inter- 
preted as "conviction" i.e., the power to retain someones 
own opinion. The second term signifies the extent to 
which somebody get influenced by another. Here both 
the conviction parameter and the influencing parameters 
are same and moreover they are identical for every indi- 
vidual. The opinion exchange for the generalised LCCC 



In this section we will discuss about the percolation 
behaviour of the geometrical clusters formed on a square 
lattice LCCC model. Here we assume that the agents 
are placed on the sites of a square lattice and follow the 
LCCC dynamics. We deflne a geometrical cluster as con- 
sisting of the adjacent sites having opinion value more 
than or equal to a prefixed threshold opinion value (f2). 
In our numerical simulation we have used random sequen- 
tial updating rule. For each value of A and fi, when the 
system reaches a steady state we measure the percolation 
order parameter Pmax = Sl/L"^ (where 5*/, is the size of 
the largest cluster and L is the linear size of the system) . 
The value of Pmax increases with A and at some A^ the 
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FIG. 2: (Color online) Comparative plots for the largest clus- 
ter size with conviction parameter for three different system 
sizes and at three various values of the opinion threshold 
{fl = 1.0,0.80 and 0.60). 



system undergoes a percolation transition (Fig. [IJ. The 
value of decreases with decrease in 51, approaching the 
value Ac as — > 0.0 (Figs. [2] and [3]). Moreover it is also 
evident from Fig. [2] that the finite size effect diminishes 
with decrease in fl. 

The percolation transition is characterised by power- 
law variation of different quantities. The order parameter 
which means the relative size (Pmax) of the largest cluster 
varies as 



(A? - A) 







(3) 



and the correlation length diverges near the percolation 
transition point as 



e^(A?-A)-^ 



(4) 



where, A^ is the critical conviction parameter. The values 
of the critical exponents (3 and v specify the universality 
class of the transition. 

However, the exponents are not determined from these 
definitions due to finite size effects. The critical expo- 
nents arc determined from the finite size scaling relations 
[l2, 30]. For example, the order parameter is expected 
to follow the scaling form 



P,r 



L"" (A? - A) 



(5) 



where is a suitable scaling function. If we plot 
PmaxL^^'^ against A for different system sizes but fixed Q, 
then by tuning P/v, all the curves can be made to cross 
at a single point. The value of A for which this happens 
is the critical conviction parameter (A^). To estimate z/, 
PmaxL^^" is to be plotted against (A^ - X)L^/'' and by 
tuning l/v, the curves are made to collapse, giving an ac- 
curate estimate of the exponent v. The other exponents 
can be obtained from scaling relations 



FIG. 3; (Color online) Plot of the critical conviction param- 
eter (X^) with the threshold opinion value (Q). 
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FIG. 4: (Color online) PmaxL^^'^ plotted against the convic- 
tion parameter A for 57 = 1.0 and fi = \. The curves for dif- 
ferent system sizes (L = 60, 100, 200, 400, 500 and 700) cross 
at = 0.760 ± 0.001 for P/u = 0.130 ± 0.005. In the inset 
the data collapse for Pmax with (X^ — A) has been shown for 
Q. = 1.0 giving l/u = 0.80 ± 0.01 and P/u = 0.130 ± 0.005. 



For il = 1.0, we plot PmaxL^^'^ against A 
(Fig. d]). The curves for different system sizes {L — 
60,100,200,400,500 and 700) cross at a point when 
iS/v = 0.130±0.005 and the crossing point (A^ = 0.760± 
0.001) gives the critical conviction parameter. Now to de- 
termine jy, we plot PmaxL^^'^ against (A^ — A)L^/'' and 
by tuning the value of l/v all the three plots are made 
to collapse on a single curve (inset of Fig. |4]) giving an 
estimate of l/iy — 0.80 ± 0.01. Although with decrease 
of f2, the critical point for percolation approaches Ac, 
the exponents remain same. We have shown the same 
plots for n = 0.80 in Fig. [5l The corresponding critical 
point is XP = 0.6955 ± 0.0005, but /3/i/ = 0.130 ± 0.005 
and = 0.80 ± 0.01. The values of P/i^ and l/i' are 
different from that obtained for the percolation tran- 
sition in case of static Ising {Ps/vg — 0.052 ± 0.002, 
l/vs = 0.996 ± 0.009) [13, dynamic Ising {(id/vd = 
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FIG. 5: (Color online) PmaxL^^'^ plotted against the convic- 
tion parameter A where f2 = 0.80 and /x = A. The curves 
for diflFerent system sizes (L = 60, 100, 200, 400, 500 and 700) 
cross at Ag = 0.6955 ± 0.0005 for filv = 0.130 ± 0.005. In the 
inset the data collapse for Pmai with AJ! — A has been shown 
for n = 0.80 giving Ijv = 0.80±0.01 and Pjv = 0.130±0.005. 
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FIG. 7: (Color online) Fourth-order reduced Binder cumu- 
lant of percolation order parameter (Pmax) for six different 
system sizes (L = 60,100,200,400,500 and 700) at ^ = 1.0 
and /X = A; the crossing point determines the critical point 
(Ac = 0.760 ± 0.001). The critical Binder cumulant value is 
IJ* = 0.62 ± 0.01. Inset shows the data collapse for the same 
value of Ijv as obtained for the data collapse of Pmax ■ 
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FIG. 6: (Color online) The cluster size distribution for LCCC 
model for three different system sizes at f2 = 1.0 and corre- 
sponding critical conviction parameter A*! = 0.760. All the 
curves decay algebraically with an exponent 1.82 ± 0.01. 



0.20 ± 0.05, = 0.85 ± 0.05) pf and standard per- 
colation {f3/v = 5/48, l/i^ = 3/4) [1] for two dimensional 
system. 

We have also studied the cluster size distribution for 
a fixed value of VI and for three different system sizes 
(L = 200,400 and 500). For 17 = 1.0, at the critical 
point (i.e. A — 0.760) all the curves decay algebraically 
as ^(5*) ^ S~'^ (where S denotes the sizes of the cluster) 
with T = 1.82±0.01 (Fig. [6]). The value of r remains same 
for other values of VL (at corresponding critical values of 
A). 

For further verification of the critical point and the 
universality class, we have studied the reduced fourth- 
order Binder cumulant of the order parameter, defined 
as [13 



FIG. 8: (Color online) Fourth-order reduced Binder cumulant 
of percolation order parameter [Pmax) for different system 
sizes (L = 60, 100, 200, 400, 500 and 700) at 17 = 0.80 and 
fi = X; the crossing point determines the critical point (AJ? = 
0.6955 ±0.0005). The critical Binder cumulant value is U* = 
0.62 ± 0.01. Inset shows the data collapse for the same value 
of l/u as obtained for the data collapse of Pmax- 



U = l 



I pi 



(6) 



where Pmax is the percolation order parameter (as de- 
fined before) and the angular brackets denote ensemble 
average. ?7 — >■ | deep inside the ordered phase and 
U ^ Q in. the disordered phase when the fiuctuation 
is Gaussian. The crossing point of the different curves 
(t/ — A) for different system sizes gives the critical point 
(AP = 0.760 ±0.001) for Q. = 1.0, which is in good agree- 
ment with the previous estimation from finite size scaling 
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FIG. 9: (Color online) PmaxL^^'^ plotted against the convic- 
tion parameter A where = 1.0 and fj, = 0.50. The curves 
for diflterent system sizes (L = 60, 100, 200, 400, 500 and 700) 
cross at A? = 0.842 ± 0.001 for jS/v = 0.130 ± 0.005. In the 
inset we have shown the data collapse for Pmax with A? — A 
for Q. = 1.0 and = 0.50 giving 1/v = 0.80 ± 0.01 and 
/3/z^ = 0.130 ±0.005. 



FIG. 10: (Color online) PmaxL^^" plotted against the con- 
viction parameter A where = 1.0 and /i = 1.0. The curves 
for different system sizes (L = 60, 100, 200, 400, 500 and 700) 
cross at A? = 0.687 ± 0.001 for P/u = 0.130 ± 0.005. In 
the inset we have shown the data collapse for Pmax with 
Ag - A for n = 1.0 and ^l = 1.0 giving l/u = 0.80 ± 0.01 
and = 0.130 ±0.005. 



of the corresponding Vl (Fig. [7]). The value of U at the 
critical point for any value of \s U* — 0.624 ± 0.002 
(see Fig. [7]). The Binder cumulant also follows the scal- 
ing form 



(7) 



where W is a suitable scaling function. The data 
collapse for i7 = 1.0 has been shown in the inset of 
Fig. [71 and the value of 1/v is 0.80 ± 0.01 which is in 
good agreement with the value of l/v obtained from 
the finite size scaling of the largest cluster size. The 
same plot has been shown in Fig. [5] for 17 = 0.80, 
which also gives the same value of l/v, which indi- 
cates that the critical exponents are independent of f2. 



IV. PERCOLATION IN GENERALISED LCCC 
MODEL 

We have also investigated the percolation transition in 
case of generalised LCCC model in which the conviction 
parameter (A) and the influencing parameter (/i) are dif- 
ferent. We have studied the percolation transition for two 
sets of parameters : Q, — 1.0, ji — 0.50 and Q, = 1.0, ^ — 
1.0. In both of the cases the plots of PmaxL^^^ with A for 
different system sizes (L = 60, 100, 200, 400, 500 and 700) 
cross at a single point for fi/v = 0.130 ± 0.005 (Figs. [9] 
and fTU)) which is same as obtained for the LCCC model. 
The critical points are different (AJ? = 0.842 ± 0.001 for 
^ = 0.50 and \p = 0.687 ± 0.001 for /i = 1.0). The value 
of 1/z^ is obtained from the finite size scaling of the largest 
cluster size (inset of Fig. 1^1 and fTU)) and the estimated val- 
ues of P/v = 0.130 ±0.005 and l/v = 0.80 ±0.01 are the 
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FIG. 11: (Color online) Fourth-order reduced Binder cumu- 
lant of percolation order parameter {P-max) for different sys- 
tem sizes (L = 60, 100, 200, 400, 500 and 700) at = 1.0 and 
/i = 0.50; the crossing point determines the critical point 
(A? = 0.842 ± 0.001) The critical Binder cumulant value is 
U* = 0.624 ± 0.002. Inset shows the data collapse for the 
same value of 1/;/ as obtained for the data collapse of Pmax- 



same as that obtained for the LCCC model and are inde- 
pendent of the value of /i, which is in contrary with the 
results obtained for the transition of the average opinion 
value, where the critical exponents change with /x. This 
implies that the percolation transition is much more ro- 
bust than the average opinion transition. The plots for 
the Binder cumulant also satisfy the crossing point and 
the critical exponents as obtained previously (Figs. [11] 
and fTS]) . The cluster size distribution also decays alge- 
braically with an exponent 1.82±0.01 for ji — 0.50 which 
is the same as that obtained for the LCCC model. 
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FIG. 12: (Color online) Fourth-order reduced Binder cumu- 
lant of percolation order parameter (Pmax) for different sys- 
tem sizes (L = 60,100,200,400,500 and 700) at ^ = 1.0 
and ^ — 1.0; the crossing point determines the critical point 
(Ac — 0.687 ± 0.001) The critical Binder cumulant value is 
U* — 0.624 ± 0.002. Inset shows the data collapse for the 
same value oi l/u as obtained for the data collapse of Pmax- 



V. DISCUSSION 

We have investigated the geometrical percolation tran- 
sition of square lattice LCCC model and have found 
the critical points and the critical exponents {P/v ~ 
0.130 ± 0.005, 1/z/ = 0.80 ± 0.01, r = 1.82 ± 0.01) char- 



acterising the transition. Although the system does not 
show any finite size effect in case of the transition of the 
average opinion, the percolation transition shows promi- 
nent finite size effect for a given threshold opinion value 
(il). The finite-size effect diminishes gradually as we de- 
crease the value of ft. The transition point also decreases 
with D, but the change is continuous. The critical ex- 
ponents are independent of the value of the threshold 
opinion value as well as the value of the conviction and 
influencing parameter which shows the robustness of this 
percolation transition in this system. The critical expo- 
nents are significantly different from those obtained in 
case of static and dynamic Ising system and standard 
percolation. These exponents suggest that this LCCC 
model belongs to a separate universality class from the 
viewpoint of percolation transition. 
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